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SUMMARY 

A number of calculations of bending-torsion wing flutter are 
made at two Alach numbers, A. 1=0 ( incompressible case ) and 
A 1=0.7, and the results are compared. The air forces em- 
ployed for the case of A£ =0.7 are based on Frazer’s recalcula- 
tion of Possio’s results, which are derived on the assumption of 
small disturbances to the main flow. For ordinary wings of 
normal density and of low bending frequency in comparison 
with torsion frequency, the compressibility correction to the 
; flutter speed appears to be of the order of a few percent; whereas, 
the correction to the flutter speed for high-density wing sections, 
such as propeller sections, and to the wing-divergence speed in 
general may be based on a rule using the (1 — factor and, 
■for Al=0.7, represents a decrease of the order of 17 percent. 

INTRODUCTION 


for values of tlie reduced frequency k much beyond 1 . Thus, 
the transonic or supersonic ranges of speeds are not con- 
sidered in the present paper. The purely supersonic case 
for small disturbances is also tractable, and flutter calcula- 
tions for this case are being prepared. 

PROCEDURE 


In the idealized case of a wing with two degrees of freedom, 
wing bending and wing torsion, and based on two-dimen- 
sional air forces, the determinant al equation yielding the 
flutter condition may be written in the form (references 4 
and 5) 


Raaf~ H aa 
\R ca~ i eft 


Pahf~Hah 

11 rh ]~ eh 


= 0 


The question of the influence of the compressible properties 
of a gas on wing flutter is, of course, directly tied to the 
primary problem of determining the air forces and moments 
on oscillating airfoils moving at high forward speeds. This 
problem has been attacked by Possio (reference 1), along 
lines indicated by Prandtl, by a procedure utilizing the 
pressure or acceleration potential and the method of lineariza- 
tion of the equation satisfied by the acceleration potential 
for small disturbances to the main flow. A review and 
summary of Possio’s work with certain simplifications have 
been given by Frazer (reference 2). Frazer and Skan 
(reference 3) listed improved numerical tables of Possio’s 
results and made some numerical applications to the flutter 
problem. 

It appeared worth while to perform additional calculations 
along similar lines utilizing the notations and parameters 
more familiar in tills country. The present report has for 
its limited objective the reporting of the results of a number 
of pertinent calculations on bending-torsion flutter for a 
stream Mach number Al= 0.7 and the comparison of these 
results with those given in reference 4 for the incompressible 
case based on the Theodorsen theory (reference 5). 

The numerical accuracy of the results of Possio’s theory 
and method deteriorates as M approaches unity and as the 
frequency increases. It has been estimated that the theory 
is not safely applicable much beyond M=0.7, nor at 2l/=0.7 


The two real equations contained in the complex determi- 
nant may yield in any given problem the two unknowns, 
flutter speed and flutter frequency. 

Expressions for the R ’ s and /’s in the incompressible 
case are listed following the definitions of the various symbols 
in the appendix, and the evaluation of the terms is facili- 
tated by the use of table I. In the compressible case the 
R’s and 7’s are expressed in terms of the notation of Frazer 
and Skan (reference 3), and table II contains values for 
A/=0.7. 

Instead of a direct solution of the flutter speed and flutter 
frequency for a case in which the wing structural parameters 
are given, it is more convenient to solve for the parameter 
oj a /w 0 (the ratio of frequency in bending to frequency in 
torsion) which belongs to the border-line case of flutter for a 
definite chosen value of the parameter l/k. The elimination 
of X from the two equations contained in the determinant 
yields a quadratic equation in («*/“<»)* from which «*/«„ may 
be found and, subsequently, X may be evaluated ; and X, 
together with the given value of l/k, determines the flutter 
speed and frequency. It is convenient first to perform the 
calculations for the compressible case with A/=0.7 and 
finally, in order to furnish the desired numerical comparisons, 
to perform the calculations for the incompressible case (M=0) 
utilizing the given structural parameters and the derived 
values of the frequency ratio u h /ca a . 
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RESULTS AND DISCUSSION 


The main numerical results axe summarized in table III 
and are shown plotted in figures 1 to 4. (No tabular values 
are included for fig. 4.) The ordinate in the figures is the 
flutter-speed coefficient vjbu a where bu a represents a con- 
venient reference speed. The abscissa is the frequency ratio 

The parameter k may be considered to determine the 
wing density at a given altitude; thus, k= 0.025 represents 
the highest wing density used and *=0.2, the smallest. 
Alternatively, the change in k may be interpreted to repre- 
sent a change in altitude for a given wing, and a change from 
ft=0.10 to k= 0.05 may be considered to represent an altitude 
change from sea level to an altitude at which p equals one- 
half the density at sea level, or approximately 18,000 feet. 
The normal drop in sonic speed with altitude should be 

taken into account in interpreting M=— (see reference 3, 

c 

fig. 2). The examples treated may be further classified by 
values of the parameter a representing positions of the tor- 
sional elastic axis; thus, a=~ 0.4, <z=— 0.2, and a = 0 repre- 
sent, respectively, elastic axes at 30 percent, 40 percent, and 
50 percent chord from the leading edge. Also « a =0.2 repre- 
sents a position of the center of gravity 10 percent of the 
ehord behind the elastic axis. 

The figures show flutter curves calculated for the two 
values (M=0 and Af=0.7) representing respectively a low- 
speed or incompressible case and a high-speed or compress- 
ible case. For the usual circumstance of low values of the 
frequency ratio u h /u a , the effect of compressibility on the 
flutter speed is seen to be relatively small ; for the lower wing 
densities the effect is a small increase, whereas for the highest 
wing density used the effect is a small detrimental one. 

For the divergence speed (frequency «-»0} the formulas 
of the static case are applicable and the slope of the lift 
curve increases according to the Glauert-Prandtl rule, which 
yields the approximate formula 


v 


»(1 -M 2 ) 1 * 



flutter approach os ; that is, the low frequency or the static 
case is approached and the (1 — M 1 ) 1U rule given in reference 
4 appears applicable. The empirical formula for the flutter 
speed of reference 4 (p. 17), which is valid for high wing density 
and low values of « A /a) a , may be modified to read 


lL 

bu a 


(1-M*) 1 * 


V* 


2 +*+*« 


Since the slope of the lift curve does not increase in accord- 
ance with the Glauert-Prandtl formula beyond a certain 
value of Af< 1, the formula is clearly inapplicable beyond a 
certain value of M. This value may be taken roughly to be 
in the order of M=0.7 to 0.75. For ilf=0.7, the formula 
indicates a decrease in the flutter speed of approximately 
17 percent. 

The effects of internal damping and of the modes of vibra- 
tion have been omitted in the calculations. Inclusion of 
these effects would tend to reduce further the differences 
between the numerical results for the compressible and 
incompressible cases. This statement is borne out by the 
results of reference 3 for a tapered wing. 

Calculations for a wing with an aileron cannot be made 
by the method of reference 1 without very extensive and 
difficult computations. It may also be remaiked that the 
numerical tables of references 1 and 3 show the need of 
additional extensions and recalculations. 


CONCLUDING REMARK 

The main conclusion to be derived from study of the 
numerical flutter calculations is that the effect of compressi- 
bility on the flutter speed (wing bending — wing torsion, no 
aileron) for subsonic speeds with no shocks, although com- 
plicated, is relatively small in the usual cases and for a 
Madu number of 0.7 can be allowed for by corrections of 
small order to the incompressible-case results. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., December 18, 1945. 


For very heavy wings (k-» 0) the values of 1 jk for 
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Piques a.— The flutter coefficient t/baa against frequency ratio ra/o. for Tarlons values 
of * for M-0 and Af-0.7. Elastic axis at 40 percent chord; r«*- 0.25; a- -0.2; r«=0.2. 


Piques 4.— The flutter coefficient against frequency ratio atja, for two values of k for 
M-0 and Af-0.7. Elastlo ails at 30 percent chord; r«*-0.25; a- -0.4; r«-0.1. 
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TABLE I 

THE R' AND V FUNCTIONS FOR M = 0 
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1/* 

a 
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!.»' 
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l.k' 



-0.6 

" 

— 0. 37500 

-0.50000 
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10 
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0 
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A 
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0 
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0 
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0 
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TABLE II 

THE R’ AND /' FUNCTIONS FOR M=0.7 
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TABLE III 

NUMERICAL RESULTS 
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1.691 

2.62 

.234 

.62 

1.62 


3.860 

0 

.460 

1.740 

2.99 

0 

.57 

1.72 


1.400 

1.000 

1.070 

1.500 






2.000 

.663 

.773 

1.647 

L79 

.653 

.86 

1.54 

.100 

2.600 

.519 

.666 

1.666 

223 

.519 

.77 

1.71 


3. S3* 

.361 

.554 

LS47 

2.85 

.361 

.67 

1.&0 


. 4.960 

0 

.430 

2.130 

260 

0 

.50 

2.16 


1.700 

1.100 

1.150 

1.950 

1.13 

1.100 

1.27 

1.44 


2.000 

.818 

.601 

1.803 

1.77 

.818 

1.00 

L77 


2.000 

2.002 

1.970 

3. Ml 

. 





2: 500 

.690 

.757 

1.893 

2. a 

.690 

.89 

1.95 


2.500 

6.106 

5-785 

14.462 

_ 





3.333 

.494 

.627 

21091 

3.1 1 

.494 

.73 

2.28 


6.000 

.292 

.490 

2. 448 

4.30 

.292 

.61 

2.63 


6.490 

0 

.410 

2.630 

5.40 

0 

.64 

2.90 


2.600 

.883 

.943 

2.368 

2.80 

.889 

1.04 

2.39 


2.500 

1.297 

1. 295 

3.237 

1.81 

1.297 

1.42 

2.66 


3.333 

.638 

.724 

*21412 

8. 34 

.638 

.81 

2.72 


3.333 

1270 

2.120 

7.063 






5-000 

.436 

.560 

2.709 

4.~8Q 

.435 

.60 

3.19 


0.860 

0 

.350 

3.420 

7.84 

0 

.51 

4.00 



a 

—0.2 






1.030 

1.000 

1.050 

1.QS0 

as3 

1.000 

L 17 

0.62 


2.000 

.446 

.622 

1.243 

L 49 

.446 

.77 

1.14 


2.000 

3.242 

3. 107 

0.215 






2.600 

.269 

.516 

1.2S8 

1.81 

.209 

.68 

L23 


2.500 

4.791 

4.545 

11.303 






2.0SO 

0 

.440 

1.310 

2.09 

0 

.si 

1.28 


1.350 

1.000 

1.030 

1.390 

.65 

L000 

1.21 

.79 


2.000 

.574 

.698 

1.396 

1.62 

.574 

.84 

1.36 


2.000 

2.365 

2.251 

4.501 





.100 

2.500 

.435 

. 596 

1.487 

1.99 

.435 

.74 

1.48 


2.500 

3.922 

3.649 

9.122 






3.333 

.253 

.480 

1.599 

2 49 

.263 

.64 

L 61 


. 4. W0 

0 

.410 

1.650 

289 

0 

.58 

1.68 


' 1.700 

1.000 

1.000 

1.690 

1.22 

1.000 

1.09 

1.8S 


2.000 

.766 

.828 

1.057 

1.87 

.755 

.96 

L60 


2.000 

1.347 

1.323 

2.647 

1.25 

1. 347 

1.44 

1.79 


2.600 

.694 

.696 

1.737 

220 

.,594 

.82 

L80 

.050 

2.600 

1.728 

1.037 

4.094 






3.333 

.431 

.567 

1.890 

286 

.481 

.70 

2.66 


3.333 

2.641 

2.408 

8.027 




__ 


5.000 

.204 

.423 

2.115 

3.83 

.204 

.68 

2.21 


5.700 

0 

.370 

2.130 

4.26 

0 

-M 

2.29 


2.500 

.838 

.876 

2.189 

2.27 

.838 

.97 

2.20 


2.500 

L050 

1.047 

2.019 

L 96 

1.050 

1.15 

2.24 


3.333 

.588 

.666 

2. 221 

3.24 

.688 

.77 

2.49 

.025 

3.333 

1.459 

1.364 

4.547 




_ 


6.000 

.308 

.503 

2.517 

4.65 

.368 

.61 

2.83 


6.000 

2.279 

2.010 

10-048 




_ 


8.430 

0 

.340 

2.870 

6.46 

6 

.49 

3.17 

«=0 


1.160 

L000 

1.020 

1.170 

0.12 

L000 

1.16 

0. 14 


1.500 

.600 

.710 

1.000 

.85 

.600 

.90 

.77 


2.000 

.343 

.543 

1.086 

1.23 

.343 

.77 

.94 


2.000 

1.406 

1.410 

2.832 





0.200 

2.500 

.120 

.444 

1.110 

1.45 

.126 

.70 

1.01 


2.500 

1.518 

1.458 

3 644 






3.333 

1.570 

1.500 

5.000 






5.000 

1.590 

1. 512 

7,560 






2.610 

0 

.430 

1.110 

1.48 

6 

.69 

1.02 


1.350 

1.000 

1.010 

L 300 






1.500 

.800 

.850 

L 270 






2.000 

.615 

.637 

1.276 

1.43 

.515 

.83 

1.18 


2.000 

1.416 

1.352 

2.703 


_ 

_ 


.100 
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1.6$ 

.700 

.90 

1.51 


2.000 

1.092 

1.070 

2.141 

1.43 

1.092 

1.21 

1.73 


2.500 

. 649 

.642 

1.606 

2.13 

.549 

.78 

1.60 


1500 

1.268 

L 204 

3.009 

2.60 

1.288 

1.15 

2.99 

.050 

3.333 

.386 

.617 

1.724 

2.74 

.385 

.66 

1.81 


3.333 

1.475 

1.301 

4. 536 






5.000 

.147 

.373 

1.867 

260 

.147 

.55 

1.96 


5.000 

1.757 

1.676 

7.879 






, 5.430 

0 

.340 

L870 

3.91 

6 

.51 

L99 


2,600 

.797 

.819 

2.047 

2.36 

■■797 

.90 

2.12 


1500 

.909 

.906 

2.264 

2.31 

.009 

.97 

2.24 


3.333 

. 618 

.618 

2.059 

3.21 

.548 

.71 

2.28 

.025 

3.333 

L 174 

1.101 

3.669 






5.000 

.343 

.468 

2.288 

4.44 

.343 

.57 

2.52 


5.000 

1.426 

1.285 

0.424 






7.940 

0 

.310 

2.450 

0.67 

0 

.42 

2.76 


APPENDIX 

SYMBOLS AND NOTATION 

b half-chord used as reference unit length 

a coordinate of axis of rotation (torsional axis) 

measured from midchord (reference 4) 
p air density 

k ratio of mass of cylinder of air of diametei equal to 

chord of wing to mass of wing, both taken for 
equal length along span (This ratio may be 

expressed as 0.24 (Jff) where IF is 

weight in pounds per foot span, & is in feet, and 
p/po is ratio of actual air density to standard 
density at sea level.) 

x a location of center of gravity of airfoil measured 
from a (reference 4) 

r a radius of gyration of airfoil referred to a (refer- 

ence 4) 

u a natural angular frequency of torsional vibrations 

around a in vacuum 

u h natural angular frequency of wing in bending in 

vacuum 

d speed of forward motion 

c speed of sound in undisturbed medium 

M Mach number (»/c) 

v r flutter or critical speed 

a) angular frequency of wing vibrations 

k reduced frequency 

F and 0 functions of 1c (see references 4 and 5) 

\=2k 

In each of the following formulas the incompressible- 
case formula is listed first in the notation of references 4 
and 5 followed by the conversion of the formula to make 
use of the notation and' numerical tables of reference 3. 
(See tables I and II.) 

3f..-G.X-^-(|+«=)+(I-a-) f -(!+«) f 

=n a X- r -f+l [if,-! (i+a) (jUi+Z,-| ^)]+p Mi 

7+»+G+“)¥ 

-- 7+i[ j/ ‘-K!+“) 2 ‘] 

9 p 

"F 

=-^+p( z --H 

X a -o^4-(i+f) 

^-+p Z x 


B c 


Xa, (\ \2Q,2F 

■ = "7 +fl Ar a ji f 
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7 ~=s [-(!+“) f-(H) 2r +l-“] 

=P [^.-4(5+“) (^>+ z ‘-| z>)]+p 

-p[^-KI + “)^] 
7 ““1 [t + G _ “) 2F+1 ] 
-p( z ‘-l z ') 

J ’>-1 * F 



The functions tabulated for convenience in tables I 
and II are the primed quantities given in the following 
definitions: 


R aa — Rtui'-h®aX—^- 


Rah — R 



■p D r 

c a 

R eA = Rch + — ~ 
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